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ABSTRACT 


The equations of planar motion for a long, flexible free-free 
beam in orbit are developed and include the effects of gravity- \ 

gradient torques and control torques resulting from point actuators £ 

located along the beam. The actuators control both the orientation 
and the shape of the beam. Two classes of theorems are applied to 
the linearized form of these equations to establish necessary and 
sufficient conditions for controllability for preselected acuuator 
configurations. It is seen that the number of actuators, if pro- 
perly located, can be less than the number of modes in the model. 

After establishing the controllability of the system, the feedback 
gains are selected: (i) based on the decoupling of the original 

coordinates and to obtain proper damping and (ii) by applying the 
linear regulator problem to the individual modal coordinates sepa- 
rately. The linear control laws obtained using both techniques 
are then evaluated by numerical integration of the nonlinear system 
equations. Numerical examples are given considering pitch and vari- 
ous number of modes with different combination of actuator numbers 
and locations. The independent modal control concept used earlier 
with a discretized model of the thin beam in orbit is reviewed for 
the case where the number of actuators is less than the number of 
modes. It is seen that although the system is controllable it is 
not stable about the nominal (local vertical) orientation when the 
control is based on modal decoupling. An alternate control law not 
based on modal decoupling ensures stability of all the modes. 
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I. Introduction 


The present grant represents a continuation of the effort attempted 
in the previous grant year (May 19/7 - May 1978) and reported in Refs. 

1 and 2.* In Ref. 1, a discretized planar model of a free-free beam 
in orbit was developed assuming the beam to be represented by a max- 
imum of three point masses connected by idealized springs which ac- 
counted for the structural restoring effects. First order effects of 
gravity-gradient torques were included. It was aspumed that two of 
the discrete masses were at the ends of the beam and that the third 
mass was at an interior position, later taken to be at the middle, when 
the beam was undeflected, and along the local vertical. Control was 
assumed to be realized by the action of one or two actuators located at 
the end masses, and implemented according to the concept of distribute 
modal control. 3 According to this concept, direct independent control 
in each mode considered is possible when the number of actuators is 
equal to the number of modes in the system model (neglecting the effects 
of higher modes not included in the model) ; when the number of actuators 
(P) is less than the number of modes (N) direct control of P modes may 
be Implemented by proper selection of control law gains and the remain 
ing N-P modes are effected according to the residual coupling in the 
control inf luence matrix according to the gains selected for the P 
actuators. 

In Ref. 2, a mathematical model of a long, flexible free-free beam 
in orbit was obtained using the formulation developed by Santini (in 
modified vector form) which develops the general equations of a flexible 
spacecraft in a gravitational field. The motion of a generic point in 
the body is described as a superposition of rigid body motion plus a 
combination of the flexible structural modes. The beam's center of 
mass was assumed to follow a circular orbit, the beam considered to 
be long and slender (shear deformation and rotational inertia effects 
neglected), and the axial deformation was assumed much smaller than the 
lateral deformation due to bending. In Kef. 2, the emphasis was placed 
on tha analysis of the uncontrolled dynamics of this system where motion 
was restricted to occur only within the orbit plane; the equations of 
motion consisted of: a single equation describing the ir plane (pitch) 

libration (rigid body rotational mode) and "n" generic modal equations 
expressed in terms of the vibrational modal amplitudes as the variables. 
For planar motion with only flexural vibrations, it was seen that the 
pitch motion was not influenced by the beam's elastic motion. 


*For references cited in this report, please see reference 
list after each section. 
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For large values of the ratio of the structural modal frequency to the 
orbital angular rate, the elastic motion and pitch were decoupled; for 
small values of thi3 ratio, the elastic motion was found to be governed 
by a Hill* s three-term equation which could be approximated by a Mathieu 
equation, and the resulting stability considered by means of a Mathieu 
stability chart. Numerical simulations verified the possibility of 
vibrational instability for very long flexible beams in near-earth 
orbits. 

In this report the control of an orbiting beam based on the 
continuum model of Ref. 2 with point actuators along the beam Is ex- 
amined. The format adapted in preparing thi3 report is as follows: 

Two papers to be presented at the following conferences respectively, 
form the bases for Chapters II and III: 

1. Second AIM Symposium on Dynamics and Control of Large 
Flexible Spacecraft, June 21-23, 1979, Blacksburg, Va. 

2. 1979 MS/ALAA Astrodynamics Conference, Provincetown, Mass., 

June 25-27, 1979 (only the contributions by A.S.S.R. Reddy 
and P.M. Bainum are Included here). 

The first paper is concerned mainly with the modelling of point 
actuators, controllability conditions for a preselected set of actua- 
tors and a sample numerical case with one actuator and pitch plus two 
modes in the model. The second paper describes two control gain se- 
lection techniques using state variable feedback. The first technique 
uses decoupling of the original linearized equations of motion as a 
criteria to select gains and the second one applies the linear regu- 
lator problem to the system equations expressed in the independent 
modal coordinates. 

In Chapter IV the independent modal control concept as apnlied 
to a discrete model of the orbiting beam developed earlier, ^ is 
reexamined according to controllability and stability considerations. 

References are given separately for Chapters II, III, and IV. 
Symbols are defined in the text when and where they are used. 

Chapter V describes the general conclusions together with recomen- 
dations for future work. 

The introductions of Chapters II and III provide further details 
of the state of the art of beam modelling with relevant references. 
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II. On the Controllability of Long Flexible Beam in Orbit. 


Abstract 

The equations of planar motion for a long, flexible free-free beam 
in orbit are developed and include the effects of gravity-gradient tor^ 
ques and control torques resulting from actuators assumed to be located 
at specific points along the beam. The control devices are used to 
control both the orientation as well as the shape of the beam. Applica- 
tion of two classes of theorems to the linearized form of these equations 
is used to establish necessary and sufficient conditions for controllabil- 
ity for different combinations of number/location of actuators with the 
number of modes contained in the mathematical model. It is seen that 
the number of actuators, if properly located, can be less than the num- 
ber of modes in the system model. A numerical example illustrates the 
controlled response to an initial perturbation in both pitch angle as 
well as beam shape. 


1. Introduction 

Large, flexible space systems have been proposed for future use in 
communications, electronic orbital-based mail systems, and as possible 
collectors of solar energy for transmittal to power stations on the 
earth's surface. Because of the inherent size and necessarily low 
weight to area ratio, the flexible parts of such systems become increas- 
ingly important and m some cases the entire system must be treated as 
being non-rigid. For meeting the requirements of these (and other) pro- 
posed missions, it will often be necessary to control both the geomet- 
rical shape as well as the orientation of the configuration. 

Previously the formulation of the dynamics of a general flexible 
body in orbit was provided by Santini.^ As a specific example, the 
equations of motion for an uncontrolled long, flexible uniform free- 
free beam in orbit were developed using a slightly modified version 
of the Santini formulation. ^ The motion of a generic point m the body 
was described as the superposition of rigid body motion plus a combina- 
tion of the elastic modes. Further it was assumed that the system center 
of mass followed a circular orbit and that the pitch (rotation) and flex- 
ural deformations occurred within the orbital plane; also the elastic 
motion was assumed to be the result of only flexural vibrations. 
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The equations were linearized about a postion of zero structural deforma- 
tion and alignment of the beam along either the local vertical or orbit 
tangent* It was seen that, in the absence of control, for small ampli- 
tude pitch, the pitch equation was uncoupled from the generic modal equa- 
tions and that the generic modal equations were dynamically coupled with 
pitch only through a second order velocity term. Numerical simulations 
verified the possibility of vibrational instability for very long flex- 
ible beams in near-earth orbits. 4 

In the present paper; the uncontrolled system considered in Ref. 4 
will be modified to include the effect of actuators located at specific 
point locations along the beam (Fig. 1). The modelling of actuator 
forces will be restricted to the case where the elastic displacements 
remain small as compared with typical beam dimensions of the order of 
hundreds of meters. For preselected sets of control devices (number and 
location) and the total number of modes in the system model, controllabi- 
lity conditions will be examinee and some representative numerical re- 
sults showing the controlled response of an initially perturbed system 
will be discussed. 


«( 


-£ 



6(x,t) = l A (t) $ n (x) 
n n 


Figure 1: Beam Configuration with First Mode Deflection and p 

Actuators. 


2. Mathematical Modelling 

A. Equations of Motion for a Thin Beam in Orbit 

The equations of motion for a thin homogeneous uniform beam whose 
center of mass is assumed to follow a circular orbit have been developed 

in Ref. 4. For the case where all rotations and transverse elas- 
tic displacements are assumed to occur within the plane of the orbit, 
and where the earth’s gravitational field is considered to be spherically 
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controllability of a long flexible beam in orbit 


symmetric, these equations can be reduced to (Eqs. (24) - 
d 2 0 2 

jp- + (3u c sin 20) /2 = T p = N p /J 

+ f “n - w c < 3si " 2 e-D - - «o c ) 2 j A n = E n /M n 


(25) of Ref. 4): 

( 1 ) 


( 2 ) 


where 


0(t) represents the pitch angle between the undeformed 
longitudinal axis and the local vertical 

A n (t) is the modal amplitude of the n^- generic mode 

u - is the n — modal natural frequency 

“c 
T 

P 

E 

n 
M 


is the orbital angular rate 
is the external pitch acceleration, N /j 

p 

is the effect of external forces on the 
generic mode 

is the generalized mass of the beam in its n^- mode 


th 


Jared S wfS^hI ^ ^ elaStlc "s?lacenents are small as com- 

the baam length- It can be concluded that there is no first 
b “ 6 ^ the elastic motion on the rigid body pitch mo .ion, 
but that the pitch motion axfects the elastic motion due to higher order 
coupling. When the ratio of structural modal frequency to orbital rate 

troneJ .r PltCh am ? Utude is —11. xt is shown" that the uncon! 
f , aid of l iltTl \° an be a PP roxlmate dbyaMathieu equation, and with 
identified. 3 U Chart P arame tric instability regions can be readily 

d , r ^ the development of the actuator modelling and subsequent con- 
sideration of controllability, Equations (1) and (2) will be linearized 
and time and length will be nondi-ensionalized according to ’ 


T = 0) t 
C 

Z n * A J Z 

n n 

where l = length of the undeformed beam. 

The resulting linearized system equations are: 

^ + 36 = T / w 2 

UT PC 

d 2 Z 


(3) 

(4; 


dx 


( Z‘ ) "n “ W w c ’ n " 1. 2, 


(5) 


( 6 ) 


i 
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By defining 


0 = X, 


and 


de/dr - x x = x^ 

dZ l /dT r X 2 " X n+2 

I 


I 


(7) 


n-1 


dZ ./dt = x„ = Xn„ 
n n-1 n in 

eqs. (5) and (6) may be written in the standard form: 

X 


where 


-A I 0 


X + B U 
c c 


(8) 


X « [ x r x 2 . . . . X n , X n+1 , . . . X 2n ] T » state vector 
0 « nxn null matrix 


I = nxn identity matrix 

r - 




o 


o 


N / Z 


B c - L 


~ 1 

L t J 

U = [T /w 2 , E /M-iw 2 , . . . 
c L p c 1 1 c 

and represents the control vector 


21 * 


E . / M - £o> 2 ] 
n-1 n-1 c J 


B. Modelling of the Point Actuators along the Beam 

It is assumed that p actuators be located along the beam at points 
g , • • gp» where 5 lies along the beam’s undetorned longitudinal 

axis and 5=0 corresponds to the mass center of the undeformed beam. 
The actual control forces associated with these actuators w^.11 be 
designated f , f 2 , . . . f j • . . fp , respectively. For small elastic 
displacements the component of the control force, f j , parallel to the 
5 axis is very small and the component parallel to the 5 axis can be 
approximated by f-i . Thus the control torque due to the j th actuator 
may be expressed by 

N = / r x f dm (9) 

f j j 


\ 
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whare 


fj " k f 6 ( 5 -^ ) 


r " 5 j i + h ( V 

? 4 - 2 ♦“ (5) A (t) k 

j n -1 2 n 


and 


( 10 ) 


<P z is the z th component of the modal shape 
function corresponding to the n th mode 
After integration there results 

Np ° -j fj£j x (const for a uniform beam) 

For convenience the constant will subsequently be incorporated into 

( 11 ) 


f y It is then clear that for p actuators. 


N 


£ ■ -3 (£ 1 5 l + f 2 «j + . . . + f p « p ] 


and that this term divided by the pitch axis moment of inertia, J, 
provides the control acceleration for the pitch motion. 

For the generic modal equations the control forces can be trans- 
formed in^o the corresponding modal forces by^>4 


E = / <}> n • f dm 

n j 3 


( 12 ) 


(13) 


Under the assumptions previously stated, 

E n = x ( ccsnst *) 

As before, the constant will be incorporated into fj so that the effect 
of all p thtusters on the n th generic mode can be expressed by 


f ! *1 V + • • • + £ p *” <y 


(14) 


The control vector, U , can now be related to the actuator forces, 
actuator locations, arid modal shape functions by 
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l * 


where T 

x " i x x » x 2 * • • • x n » x n+1 » • • • x 2n l 

I - nxn identity matrix 

B » B nxp matrix 
act 

f . [f r f 2 . . . f p ] T 

The system represented by Eq. (17) is controllable if the pair 
[A,B] is controllable 7 - i.e. it Cc'n be proven that the controllabi- 
lity matrix associated with the original state and control coefficient 
matrices 


has rank 2n if and only if the controllability matrix associated with 
the pair of reduced state and control matrices, 

[A , B) 

has rank n.^ 



Furthermore, if the matrix A has eigenvalues of unit multiplicity 
(i.e, non-repeated eigenvalues), the system given by Eq. (17) is con- 
trollable if and only if each row of B has a non-zero entry. 7 

For the case where A has repeated eigenvalues (multiplicity greater 

than one) , the reduced order state matrix. A, can be written in Jordan 

block matrix form. If the eigenvalues A^, A 2 . . . A m have multiplicity 

of ni, no, . . . n , respectively, with n.. ® n, then A canbetrans- 

_ x , *• m 1=1 

formed as 
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The system (17) can be divided into m subsystems. For the system (17) 
to be controllable, these m subsystems -ith tneir corresponding blocks 
in the B matrix must each be separately controllable. 

B. Application of the Controllability Theorems 

The theorems briefly outlined here will now be applied to several 
cases of interest for different combinations of numbers and locations 
of the actuators along the beam. 

Case 1: One actuator at one end of the beam with pitch and two generic 

modes contained in the mathematical model. 

The actuator is assumed to be located at the left end (Fig. 1) £ * 
-i/2. The first and second modal shape functions for a free-free beam 
can be evaluated at that point to yield® 


*1 (- J ./ 2 ) - $1 (- 1 / 2 ) = 2 


and Mi - M 2 - 

The state equation for this system can be expanded in the form of Eq. 
(17) with the result that 


|*3 0 0 

j 0 c 0 

0 0 d 


where 


i/2Jco 


c = (ojj/u^)' 


b » 2/M^w* 
d = (oJ 2 /w c )' 


The controllability matrix based on the reduced system matrices, A 
and B, becomes 


C = [B ! AB ! A 2 B] = 


3a 9a 
be be 2 
bd bd 2 


For controllability the matrix, C, must have a rank of or 
det C = - ab 2 (c-d) (c-3) (d-3) i 0. 

The necessary and sufficient conditions for controllability become 

i u) 0 (trivial) ; u. t /3u> ; w. / /3w 

i z ± c z c 

The last two of these conditions will, in practice, place a lower bounds 
on the stiffness and/or an upper bounds on the length of such a long flex- 
ible structure in orbit. 
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Case 2: Three actuators two of which are assumed to be placed at 

the ends and one at the mid point of the beau, with the 
mathematical model containing only the first two generic modes. 

For this case. 


c 0 
0 d 


a ■ b * *^2/M^iw^ with - M 2 

and c,d defined as in Case 1. 

The controllability matrix may be calculated as 

I 

^ a a -b | ac ac -be 
_a -a 0 j ad -ad 0 

It can be seen that since the B matrix itself has rank 2, then C 
will automatically have rank 2 and the system concrollability is inde- 
pendent of the nature of the matrix A. 

Case 3: Two actuators one each at the ends with pitch plus the first 

generic mode in the system model. 

For this system the A and B matrices in Eq. (17) become 

r- r- —1 


1/ 2Jw^ 
c 


a -a 
b b 


2/Mj^ 


(w./td )‘ 

1 c 


The resulting controllability matrix 

„ a -a I 3a -3a 
C - I 

b b i be be 

has rank 2 since the B matrix has rank 2. Thus, the system controllabi- 
lity is ensured. 

Case 4: Three actuators two of which are assumed to be located at 

the ends and the remaining one at the mid point of the beam; 
the model contains pitch plus the first two generic modes. 


For this case 


j*3 0 0 

I 0 c 0 

! 0 0 d 


-a 0 


-b 0 
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with c,d as defined in Case 1 
and a « £/2Ju 2 ; b » 2 


and a « £/2Ju 2 
and H « M 0 ° 


2/M, Aw 2 
i c 


/2/M, Jc.o 2 
i c 


* 


a -a 0 

1 3a 

-3a 

0 1 

9a -9a 

0 


c - 

b b -e 

! * 

be 

-ce 1 

1 

be 2 be. 2 

-c 2 e 



b -b 0 

| bd 

-bd 

1 

0 ) 

bd 2 -bd 2 

0 


must have 

i rank 3 to 

ensure 

the 

system 

is control 1 able, which means that 


from the nine columns it must be shown that at least one 3x3 aon-zero 
determinant exists. 

If we arbitrarily select the first, third, and fourth columns, 
then -eab (d-3) t 0 

which is guaranteed if C 02 ^ ^3o) c . Although this is a sufficient con- 
dition for controllability at this point we don't know whether it is 
also a necessary condition. 

As an alternate, let us select the first, third, and sixth columns 
of Cj then 

ab 2 (d-c) (c-3) (d-3) 0 

which will be ensured if the following sufficiency 
conditions are satisfied: 

w 2 ^ w i* f ^3o » and 

In order to establish necessary conditions, we will now assume that 
any two of the frequencies are the same and then apply the theorem for 
the case of repeated eigenvalues. 

(a) First if it is assumed that c = d =* u^), the corresponding 
subsystem matrices are 


'c 0 
0 c 


b b 
b -b 


Since the B matrix has rank 2, the condition: c ^ d, is not 

necessary for controllability. 

(b) If it is assumed c » 3, (io^ - /3o) c ) the corresponding sub- 
system matrices are 

A- [ 3 °] B- \ a “ a °1 

LO 3_ Lb b -ej 

Since B has rank 2, the system is controllable for this case. 
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(c) We now consider the case when d ■ 3 (u- ** ) where the 

Qll^CVChom m o -r- -f o * C 


subsystem matrices are 

A " [o 3] B ’ 


a -a 
b -b 


C » a -a 0 j 3a -3a 0 j 

. b -b 0 | 3b -3b OJ 

The controllability matrix has only one independent column and 
can not have rank 3 when 

Finally we consider the case where 3 « c =d or * *■ u-j = w* 

(admittedly, of only academic interest) C1 1 




a 

-a 

0 


B = 


b 

b 

-e 




b 

-b 

0 


3a 

- 

•3a 

0 

J 9a 

i 

-9a 

1 3b 


3b 

-3e 

i 9b 

9b 

3b 

- 

•3b 

0 

I 9b 

-9b 

:hat 

the C 

matrix can 

not 

say 

that 

only 

one of 

the 


9 

0 

-9e 

0 


w 2 £ /iw c 

In actual practice the repeated frequencies associated with the modes 
^idu^ad in the model will be known so that it is possible to verify 
in advance of the selection of the control law whether the particular 
choice ^f number and location of actuators will result in a con- 
trollable system. 

4. Numerical Example 

A numerical example of Case 1 is considered where it is assumed 
that the control force generated by the single actuator depends on 
only rate feed-back according to 

f i = K i*i + k 2 x 2 + k 3 x 3 

where the xi ara the pitch, and non-dimensionalized first, and second 
modal amplitude derivatives, respectively, with respect to the orbital 
time, t. 
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free bean is^l/ioo natural fre 9 ue n c Y of the free- 

tude circular orbit. For this case C '"' m ° VeS * 250 “* ^ alti “ 
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— ° f ««*« 

of 1.06cm, oould cxhlbulhSc J'c^eccies"" 

panded h to C yleld' :Sl> ’ t>0 “ trlvIi>l P“f of Eqs. (16) or (17) may be ex- 
*1 + 3 Xi~ 59 -52K 1 x 1 -59.52K 2 x 2 -59.52K 3 x 3 “ 0 
X 2 + 3200x 2 -20.0K 2 x 2 -20. OK^ -20.0K 3 x 3 ■ 0 
X 3 +28,800x 3 -20.0K 3 x 3 -20.01^ -20. OK^ - 0 

-O.OlGgWnotfiS SeleC - K 1 = “°- 00577 * x 2 " -°«05656, and K, - 
j . . ^ ese E a ms would correspond to much less than critical 

?hTSe“ .PPeart aoTao^e 

iolflal X (0, - 0, Cbe L'iSSUi^U,'^ 

2L*ss r’L u-*T2sr tlne " ith ths reiativeiy io " ^ 



Figure 2: Case 1 - Controlled Response, Pitch + Two Modes 

with One Actuator at Left End. 
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Figure 2: Case 1 - Controlled Response, Pitch + Two Modes 

with One Actuator at Left End. 


This example is presented as a verification that the system in 
Case 1 with the number of actuators less than the number of modes is 
corcrollable. In a related paper methods of selecting control law gains 
based on decoupling considerations is discussed. 9 Control gains are 
selected based on the following two criteria: (i) decoupling of the 

linearized system equations with appropriate state variable feedback; 
and (ii) applying the linear regulator problem to the n modal coordi- 
nates separately and thus selecting the gams by solving groups of n 
two dimensional matrix Riccati equations. y 

5 • Concluding Remarks 


In the present paper a model is developed for predicting the dyna- 
mics of a long, flexible free-free beam in orbit under the influence of 
control devices which are considered to act at specific points along the 
beam. Application of two classes of theorems establishes the necessary 
and sufficient conditions for controllability and clearly demonstrates 
that the number of actuators, if properly located, can be less than the 
number of modes in the system model* 

The insight gained by this preliminary study will be useful in 
analyzing the dynamics and control of more complicated structures such 
as of a large flexible plate in orbit, which more adequately represents 
a large flexible orbiting platform. Another possible extension to the 
current work would be a study of the effect of using control devices 
which are distributed along the beam instead of being treated as point 
actuators. 
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III. Decoupling Control oi a Long Flexible Beam in Orbit v 


Control of large flexible systems using state 
variable feedback is presented with a long flexible 
beam in orbit as an example. Once the controllabil- 
ity of the system is established, the feedback gains 
are selected: (i) based on the decoupling of the 
original coordinates and to obtain proper damping 
and (ii) by applying the linear regulator problem to 
the individual modal coordinates separately. The 
linear control laws obtained using both techniques 
are then evaluated by numerical integration of the 
non-linear system equations. The response of the 
state together with resulting beam deflection and 
actuator force (s) required are obtained as functions 
of time for different combinations of the number/lo- 
cation of actuators and the number of modes In the 
system model. Also included are results showing the 
effects (control spillover) on the uncontrolled modes 
when the number of controllers is less than the num- 
ber of modes, and the effects of inaccurate knowlege 
of the control influence coefficients which lead to 
errors in the calculated feedback gains. 

1. Introduction 


Future proposed space missions would involve large, inherently 
flexible systems for use in communications, as collectors of solar 
energy, and in electronic, orbital-based mail systems. ’ For the 
first time the flexible parts, and in some cases the entire system, 
due to its size, must be modelled as being completely flexible. In 
order to satisfy the requirements of such missions, it will be necessary 
to control not only the orientation of the system but also the geo- 
metrical shape of the configuration. 


As a specific example of the general formulation of ^ .a dynamics of 
an arbitrary flexible body in orbit developed by Santini-, the uncon- 
trolled motion of a long, flexible beam was investigated. 4 The motion 
of a generic point in the body was described as the superposition of 
rigid body motion plus a combination of the elastic modes. 
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Further it was assumed that the system center of mass followed a cir- 
cular orbit tnd that the pitch (rotation) and flexural deformations 
occured within the orbital plane. For this planar motion, it was seen 
that the pitch motion was not influenced by the beam’s elastic motion. 
The decoupling of pitch and the elastic modes was observed for large 
values of the ratio of the structural modal frequency to the orbital 
angular rate. When the values of this ratio are small the elastic mo- 
tion is governed by a Hill’s three- term equation which could be approxi- 
mated by a Mathieu equation, and the resulting stability was considered 
by means of a Mathieu stability chart. Numerical simulations verified 
the possibility of vibrational instability for a very long uncontrolled 
flexible beam in near-earth orbits.^ 

The controllability of a long flexible beam with point actuators 
located along the beam is considered in Ref. 5 for the case of small 
amplitude flexural deformations (Fig. 1). Necessary and sufficient 
conditions for controllability with preselected locations of actuators 
are derived using theorems developed in Ref. 6. Once controllability 
is assured, values for the gains in the control laws are selected on an 
arbitrary basis, and only for one combination of actuator location and 
number of flexural modes. 5 


2, C 



6(x,t) = Z A (t) $ n (x) 
n n 


Fig. 1. Beam Configuration with First Mode Deflection and 
p Actuators. 
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In the present paper selection of cortrol gains for any large 
flexible system using the following two criteria is discussed: (i) 

decoupling of the linearized system equations with appropriate state 
variable feedback 7 ; and (ii) applying the linear regulator problem to 
the modal coordinates (n) separately and, thus, selecting gains by 
solving groups of "n" two by two matrix Riccati equations. 8 *** » ii 

A long flexible beam in orbit is taken as an example with the 
model developed in Refs. 4 and 5. Gains are selected using the two 
techniques and numerical simulation of the non linear equations is 
employed to predict the responses for sample cases. The deflection 
of the controlled beam at various instants of time is also illustrated. 


2. Decoupling by State Variable Feedback 


After appropriate linearization the dynamic model for any flexible 
system can be represented by 

AX BX + CX = DU (!) 

where A is dn jum non singular matrix 

B,C are nxn matrices 
D is an nxm matrix 

X is an nxl state vector representing deflections 
in addition to the rigid body rotations. 

U is an mxl control vector 

Equation (1) can be written in more standard state space form oy defining 


X - x v X = x 2 s X L as 


J— 1 


— - 


- 


— — 

*i 


0 I 




0 


- 




4- 


*2 


-A _1 C -A _1 B 


1 '2 


A -1 D 

1— -J 


— - 

J 

— * 


( 2 ) 
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Equation Cl) takes into account the modelling of any structural damping 
inherent] v present in ihe system* Controllability of the systems repre- 
sented by Equation (2) can not be obtained using theorems developed in 
Ref. 6, unless B 0. In cases where B */ 0, the controllability matrix 
of the pair 


-A _1 C -A _1 B 


must have rank 2n. 


After choosing a state variable feedback control law of the form 
U - KX + LX 


K « mxn rate feedback gain mafrix 
L ** mxn position feedback gain matrix 
Equation (1) can be rewritten as: 

X + (A -1 B-A _3 TK)X + (A _1 C-A -1 DL)X = 0 

T 

For decoupling of the states X = (x^> X 2 » •••» X n ) 

the matrices (A“1 b-A”1dK) and (A~1 oA“1du) must be diagonal. 


i.e. A B 


~ 1 B-A _1 T)K 


A _1 C-A _1 DL » u 


where 


0 and oj 


Redefining 


A B » E , A C - F 
A _1 D = G 


we have 


F — GL 73 co 


column of the K matrix 


column of the L matrix 



4 

1 1 # 


E . f ® column of (E-t) matrix 

F^' *> i column of (F-co) matrix 

Equacions (10) and (11) can be uritcen as 2n sets of algebraic 
equations of the form: 


GK. 


GL. 


V 

T» I 

‘ i 


i « 1,2, . . n 


( 12 ) 


Consider one of the above sets of linear algebra 4 *! equations for the 
case where i « 1, 


GK X » (13) 

There are n equations and m unknowns (the exements in the first column 

of the K matrix) . The fundamental theorem for a set of n linear equations 

with m unknowrs is now applied* 2 

For a unique solution : 

Case 1: If n > m (more equations and less unknowns) the rank of G 

and the augumented matrix [G ! E^*] must be n. 

Case 2: If n ■ m (number of equations a number of unknowns) the 

rank of G and the augumented matrix [G J E £] must, be m 
(or n) - i.e. G must be non-singular. 

Case 3: If n < m (less equations, more unknowns) no unique solution 

exists. 


For non-trivial solution : 

Case 1: If n > m (more equations and less unknowns), the rank of G 

and the augumented matrix [Gj E^] £ m. 

Cate 2: If n 8 in (number of equations = number of unknowns), the 

rank of G and the augumented matrix [G [ E^ 1 ] £ m. If the 
rank = m = n, a unique solution exists. 

Case 3: 7f n < n (less equations and more unknowns) the rank of G and 

the agumented matrix [G { E^ f ] £ n. 

These conditions must be satisfied by the 2n subsystems given by 
(12) for decoupling to be implemented and, thus, dictate the choice 
of the actuators. 

For cases where the number of actuators are more or less than the 
number of original coordinates in the sy°tem, the controllability con- 
ditions and the conditions to be satisfied for decoupling may pose 
numerical (computational) problems, especially when the order of the 
system ie large. 
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When the number of actuators equals the number of original coordinates 
all controllability and decoupling conditions depend on the non 
singularity of matrix, G. This matrix is an nxn matrix and can be 
made non-singular by properly selecting the location of the actuators* 

If there are 6 original coordinates and 3 actuators, then the con- 
trollability matrix has 36 columns out of which 12 columns must be 
independent to have rank 12. The maximum number of determinants to be 
evaluated are 


36 r m 36! 

C 12 12124! 


= 1.2516775x10/ determinants. 


Assuming a L sec. computational time required for the evaluation of 
each twelth order determinant, the examination of all possible com- 
binations vould involve 347,688 hrs. of computer time. 

Specific cases, where controllability of this system is examined 
when the number of actuators differs from the number of modes in the 
system modAl are presented in Ref. 5, but only for a low order system 
(n - 3). 


3. Linear Regulator Problem 


Using modal analysis 8,9,10,11 dynamical systems represented by 


where 


Mq + Gq + Cq + Kq « u(t) 

q « n dimensional vector describing angular and elastic 
displacnents 

M,K = positive definite mass and stiffness matrice 
G a gyroscopic antisymmetric matrix 
C =* pervasive damping, either positive definite or 
semidefinite matrix 
u(t)=* control vector 


or by (b): 


Mq + Gq + Kq = u(t) 

M ® M T , G - -G T and K » K T 


can be transformed to 

(a) X + DX + FX =* Eu - u 1 


(b) X + FX = Eu « u 1 (17) 

with D and F representing the diagonal transformed matrices. 

The left hand sides of Eqs. (14) and (15) represent either 

damped or undamped harmonic oscillator These oscillators can be con- 
trolled optimally and independently with one control force for every 
independent modal coordinate. 
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The actual control, u(t), in the original coordinates, can be 
calculated from the control, u'(t), in the decoupled (modal) coordinates, 
X ± > by 

u » E ^u' (18) 

The component of the control vector, u', can be calculated as 
follows : 

The i £ k independent modal coordinate is governed by 


x i + D i X i + F i X i " U ’i 

where p „ <ji a g 0na i element of the D matrix 
F i ■* i th diagonal element of the F matrix 
Equation (19) can be written in state space form by defining 


(19) 


l 1 1 x 1 2 x L 


( 20 ) 


similar to the standard form 

X ■ Ax + Bu (21) 

tH 

A performance index f^r the i modal coordinate is defined as 

oo 2 

J i - / (x^Q^ + (u^) R i ) dt (22) 



where 


and R^ is a scalar. 


The control vector is given by 

T 

u = -R^bTs.X, where X 4 = Cx, X ± ] 
i i i x i x 12 

Where S is the symmetric matrix solution of the two dimensional 
Riccati^equation : 

-Vl - A I S i + s i b i r " 1 b I s i - 5 ■ 0 


(23) 


(24) 
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— ^ ^ - ***** 


- F i - D i 


B* - [ 0 , 1 ] 


R - R ± and Q ■» 


of solved S'SoSd'SrTSnh'ST'r^.S; W ’ ' he 


*12 ■ R i r~ F i ± ^W+ Q./R ] 

1 1 


R j r~ D i ± + 1_ (q. ? 2s r ] 

p x o i-t o 


F i s i + D i s i + JL s, s, 

22 12 *12 i 22 


defini^mSix.^ radicals are selected such that S jL is a positive 


V “ “ t S i ♦ S ] 
12 *22 


4. Numerical Examp le 

strate *" “J 11 ls "“«««> t. de«on- 

<Flg. 1) including poi« actuat'ots"^“«l ! arllet ' ”« ~del 

based on Che assumption that all rntart ke " f J°j < e£s ‘ 4 and 5 and is 

within the orbital “laL rae e c°atLns S „r 4a£ ° raatl °" a °™T ™ly 

i- “c. me equations of motion are given by 4 ’ 5 

d 2 0 T P 

•jzz + 30 - -X 

dt Jw 2 (2q\ 


d \ ( “ n 1 
d?- + ( — 

\ c/ 


z n = _n 

M £u)* 
n c 


n “ 1 , 2 , ... 


0 - pitch angle relative to local vertical 
T = a) c t, normalized time 

Z n “ A r/*» non-dimensional modal amplitudes 

1 " length of the beam 
a> c = orbital angular rate 
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After defining 




Case 1: pitch + 2 modes are considered in the model with actuators 

located at (-V2* V 2> V 4) 


— — 


— — 


— 

-3.0 • 


59.52 -59.52 29.76 


f l 

-3200.0 

B u =■ 

20.0 20.0 - 2.0 


f. 


c c 



2 

-28800.0 


20.0 -20.0 9.3 


**3 

— 






The gains are selected for decoupling anu critical damping in the 
decoupled mo*des (only rats feedbackis considered here since the 
uncontrolled system is already decoupled) . The required control 
forces are given by: 




-0.011 

-2.835 

-5.2757 

0.017 

-2.835 

3 5171 

0.0602 

0.0 

-17.5855 


The time response of pitch, nondimensionalized modal amplitudes, 
and forces required are plotted in the following figures, (Fig. 2). 

As the mode number (and frequency) increase the decay time decreases. 
The pitch takes a relatively long tins to decay since its natural 
frequency is very low and only rate feedback is considered here. The 
maximum amplitude of the forces are of the order of newtons for this 
application. 



Fig.* 2 a 

Fig. 2. Decoupled Controlled Response - Pitch + Two Modes 
with Three Actuators (~^/ 2 * ~/2> 0 - 
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(Non-dimensionulized 











Pig. 2. 


Decoupled Controlled Response - Pitch + Two Modes 
with Three Actuators ( - A /2> A /2> 


Case 2: pitch +> 4 inodes considered with actuators located at 

(-*/ 2,- A /4, 0, */4, l /2) 


A =* 


and 


B u » 
c c 


-3.0 0 0 0 0 

0 -3200.0 0 o 0 

0 0 -28800.0 0 0 


0 

0 

0 

-93079 

.50 

0 


0 

0 

0 

0 

-255331.40 





— 


- 

59.52 

29.76 

0.0 

-29.76 

-59.52 


f, 

20.0 

-1.98 - 

•14.0 

-1.99 

20.0 


f 

20.0 

-11.698 

0.0 

11.635 

-20.0 


c 

20.0 

-12.437 

14.18 

-12.37 

20.0 


I 

20.0 

-5.1266 

' 0.0 

5.01 

-20.0 


f 


f 
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For case 2, the gains are selected using both decoupling and an 
application of the linear regulator proolem to the independent modal 
coordinates. 


tfsing decoupling (pitch and the first four modes are critically 
damped) , the control forces are given by 



-1.199 

-1.423 

-789.03 

-7.578 

3365.27 

-3.329 

0.0 

-2183.39 

0.0 

9392.17 

-2.443 

4.015 

-1620.70 

-21.649 

6946.89 

-3.278 

0.0 

-2200.36 

0.0 

9397.81 

-1.669 

-1.423 

-780.55 

-7.578 

3362.45 



The response to an assumed perturbation of 0.01 in all Lne posi- 
tion coordinates, including pitch, is shown in Fig. 3. The maximum 
amplitudes of forces are of the order of hundreds of newtons. Fig. 3, 
also illustrates how the initially deformed beam is straightened out 
under the influence of the controllers. It is seen that after 36 secs 
the beam is essentially straight, but continues to exhibit a pitch 
displacement until about 4000 secs. 



- 0.01 

Fig. 3 A 

Fig. 3? Decoupled Controlled Response - Pitch + Four Modes 
with Five Actuators (-V 2> - Z /4> °» l /l>- 


> < 


* 


\ 
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(Non-dimensionalizcd) 









Normalized Deflection 



4 


The V near regulator problem is applied tc the system in the decoupled 
(modal) coordinates, where the non-dimensionalized rates are penalized 
by a factor of frequency squared as compared with the non-dimen— 
sionalized position coordinates. After solution of the five two 
dimensional matrix Riccati equations, the actual control, forces are 
given in terms of the following gain matrix. 


-0.056 

1.233*10"* o.o 

0-0 

0.0 

-0.6312 

-0.712 

-334.32 

-4.343 

1532.5 

-0.1559- 

0,0 

0.0 

0.0 

0.0 

-1.7527 

0.0 

-109 .17 

0.0 

4696.0 

-0.1144. 

3.545*10"* 0.0 

0.0 

0.0 

-1.2359 

2.007 

-310-35 

-6.021 

3473.4 

-0.1533 

0.0 

0.0 

0.0 

0.0 

-1.7255 

0.0 

-15.CO.13 

0.0 

4693.9 

-0.0545 

0.0 

0.0 

0.0 

0.0 

-0.5142 

-0.712 

390.27 

-4.343 

1631.2 





- 0.01 


Fig. 4 A 


Fig. 4. Linear Regulator Apolication Controlled Response - 
Pitch + Four Modes with Five Actuators ( -V? > “^4, 

0, */4> V 2 ). 
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The response to an assumed initial perturbation of 0.01 in all the 
position coordinates is shown in Fig. 4. As the rates are penalized 
heavily when compared to the positions, the controlled system fre • 
quencie 3 are not changed appreciably and the damping obtained in the 
individual modes is less than critical. The very small numbers and 
zeros (which are in reality < 10 _7 in the second to fifth columns of 
the position feedback portion of the gain matrix) are due to unit 
weighting of the positions in the matrix. It can be sho'm that 

the forces required have a maximum amplitude of the order of thousands 
of newtons. When compared with Fig. 3, the maximum amplitude of the 
forces reouired here are approximately two orders of magnitude larger. 
This can be explained by the fact that the model used here Includes the 
third and fourth higher frequency modes, and it has been assumed that 
all four modes and pitch were initially excited equally. 

5. Conclusion 

A technique for selecting control system gains based on the de- 
coupling of the original linear system equations of motion is pre- 
sented. This avoids use of modal analysis and does not require sys- 
tem matrices to be symmetric or skew symmetric. When the number of 
actuators is equal to the number of modes, a unique solution for the 
control gains depends on the non-singularity of a matrix ba=ed on 
(modal shape functions evaluated at) actuator locations. When the 
number of actuators is less than the number of modes and the order of 
the system is high, implementation of decoupling control may be limited 
by the computational capacity. 

The linear regulator problem can be applied to the decoupled modal 
coordinates only when the number of actuators is equal to the number 
of modes. Otherwise instead of solving n second order matrix Riccati 
equations, a 2nx2n matrix Riccati equation has to be solved. 
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IV. Comparison with Results Based on Independent Modal Control 


In Refs. 1 and 2 the control of the planar motion of a long flexi- 
ble beam in orbit was studied based on the concept of distributed modal 
control. The control forces generated based on this concept provide a 
means of controlling each system mode independently of all other modes 
as long as the number of modes in the system mathematical model is the 
same as the number of actuators. For the case where the number of modes 
(N) is greater than the number of actuators (P) independent control of 
P modes is possible, and the response of the remaining (N-P) modes 
depends on the residual coupling due to the P actuators. 


The mathematical model used in Refs, land 2 is based on a three-mass 
discretization of the jLree-free beam, with two of the masses assumed 
to be at the ends of the beam and the third mass at an interior point 
which later was selected *t the center of the undeformed beam. The 
beam was represented by two hinged cantilever-type members consisting 
of the end mass connected by (assumed) massless springs which were 
responsible for the structural restoring forces (Fig. 1). One of 
the results from Refs. 1 and 2 indicates that the beam, represented by 
two degrees of freedom, and containing a single actuator at one end, 
when given an initial perturbation, will not return to the equilibrium 
position when the control is based on the concept of independent modal 
control. In the present study (Chapter II), it is clearly shown that 
a beam with a single actuator at one end and with pitch and two generic 
modes in the model can be controlled and will return to a desired unde- 
formed alignment with the local vertical. 

In an effort to resolve this apparent ambiguity, we will return to 
the previously developed discretized model and examine both the con- 
trollability and stability of the system when F * 1 < N =» 2. 


The linearized equations of motion are [Eq. (3.18) of Ref. 1 or 
Eq. (38) of Ref. 2] : 


a b v. 

.. **■ 

b a v 2 


r c o 


a 13 M* (1 + nig); b = M*; c - 3cjg M* (2 + m^) + k 
M* =* m^/M 

°0 = m o /m 

m * mass of each end mass 
nig = mass of interior mass 

k = elastic restoring constant (=3EI/x, 3 for assumed 
cantilever members) 






I ■ length of each member (one-naif the undeforned oeam 
length) 


u>Q = orbital angular velocity 

v. , = linear deflection of each end mass 

1, z 


control forces due to actuators 


' 1.2 


Eq. (2) can be rewritten as 


, 

< S 

1 

1 

ac 

-be 


i — 

rH 

> 

1 

+ 1 .fa 

-b ' 


J 

1 

CM 

i> 

1 

F-b 2 

-be 

ac 


_ V 2. 

a*- b z , , 

L 

a 


F 

L V 2 J 


( 2 ) 


If 


'1 *1 :i A 3 A 1 

v„ - Xo » Xa “ 


'2 “ A 2 '2 " A 4 A 2 

thei\ Eq. (2) can be written in standard state space form as: 


“ -< 


- 



* 






-n 

^1 


0 

0 

1 

0 


x i 


0 

0 



0 

0 

0 

1 


x 2 

1 

0 

0 



-ac 

be 

n 

0 






x 3 


a z -b^ 

a z -b z 

u 


x 3 


a 

-b 

• 


be 

-ac 


o 


v 


-h 

a 

X 4 


F^F 

F^F 

U 

u 


x 4 


u 

a 

— 


— 



- i 




- 


Eq. (3) is now written in the form: 

0 


0 I 
-A 0 


X + 


B 


2 J 


(3) 


( 4 ) 


so that according to the controllability theorem, the system represented 
by Eq. (4) is controllable if and only if the controllability matrix, C, 
associated with the pair of reduced state and control matrices, [A,B], 
is controllable. _ In this case 

-c(a 2 +b 2 ) 
a^-b~ 


a z -b 2 


I 


a -b 
-b a 


2bc 




2abc 

F+b 2- 

-c(a 2 +b 2 ) 
a 2 -b 2 


(5) 


IV- 3 


c 


H II. J 



I 

?4 


Since Jet a J ^ 0, in general, 

C has rank 2 and the system (2) is controllable. 

If only one actuator is assumed to be present (i.e. F V2 * 0), then 
Eq. (3) can be written 


1 0 


-ac be 

a z -V a 2 -t 


0 1 


0 0 


be -ac 

a^ -b 2 a? — b^ 0 0 


The reduced order controllability matrix is: 


C *=| a/ (a 2 -^) 


-b/ (a^-b 2 ) 


and its determinant. 


-c(a?+b 2 )/(a 2 -b 2 ) 2 

2abc/(a 2 -b 2 ) 2 


det C = bc/Ca 2 -^ 2 ) 2 £ 0, since in general, a ^ b. 

Thus, the system is controllable with a single actuator present. 

The stability of system (6) will now be examined using the parti- 
cular control law used in Refs. 1 and 2. The linear equations of motion 
(1) or (2) can be transformed into the modal coordinates, q- a .d q^, 
and for the case where only one actuator is present have the form: 

q i + x i q i = u i ( ' 

q 2 + * 2 q 2 =[( a +b)/( a -b)j Ul = g (i 


= c/(a+b) 


and v x = q x +q 2 


<2 = c/(a-b) 


q r q 2 


Following Refs. 1 and 2, in accordance with the concept of inde- 
pendent modal control, the control in the modal coordinates was 
selected as 

U 1 = " f l q l “ f 2 q l 


Then Eqs* (7) and (8) can be written as: 

q l + f 2 q l + ^ " 0 (10) 

8(f 2 q l +f l q l ) + q 2 + A 2 q 2 " 0 (11) 

ai^ C cJn a h te ^ iSt ? C equation for the system described by Eqs. (10) and 
Ul) can be developed with the result: J 

(s 2+X 2 ) (s^s+f^)- 0 (12 

f^ O HK da ^ Ped . m0de / emalnS 3t fre< l uenc y and is not affected by the 

turba^on 8 ^ 1 ^’ ^ ^ f 2’- the contro1 removes the initial per- 

h,Ic 1 // q l’ " Senerax, the system will continue to oscillate at 

q m0 . dal frequency. The system is unstable about 

xaw fo? ft,° (strong) sense of Routh-Hurwitz where the control 

xaw for the single actuator has the form of Eq. (9). An exauole is 
illustrated by Fig. 9 of Ref. 2 for this case where f - f'lO.nd 
demonstrates the basic phenomenon of control spillover. 2 * ’ 

Instead of selecting the control law based on the independent -,i 

Sployed'^'tSTS:' 1 ' 8 ' * C “ Pled ““ f " db “ Ck la “ *• 

“l "Vl * *2<2 (13) 

Eqs. (7) and (8) can then be expressed as: 

’l + K l’l + A l’l + >V' 2 " » (14) 

3 K l’l + 1 2 + 8K 2 4 2 + 1 2 , 2 - 0 (15) 

with the associated characteristic equation 

S 4 + + (X 1+ A 2 ) S 2 + (X 2 K 1+ X l8 K 2 )s + X 1 X 2 - o (16) 

If the rate feedback gains, K and K_, are selected to be positive 
Hurwitz c^iteria^noMn^thatl ^a^/U^T" 8 ^ ^ “ 

Numerical Example 

equall^divid^w' 1 3n ? 2 ’ the t0tal maSS is selected to be 1000 kg, 
Thus^ 17 d d d between the two end masses and central mass, m Q . 8 

M = 1000 kg 
m = M/3 =* 

n o 3 m (/ m = 1 
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i + 


K* » m 2 /M * 111.11 kg 
a ■ K*(14 h1q) =» 222.22 kg 
b - M* » 111.11 kg 

k =* 3EI/Z 3 = 0.18497 N/m for a cylindrical 

wrought aluminum tubular beam, 100m. 
long (L = 2£) 

w 0 « 1.115 x 10” 3 rad/sec 

c * 0.1862522 

X 1 - 5.5876X10" 4 

X 2 *> 1.67627xl0“ 3 

g ■ 3 

If the feedback rate gains are selected as 


K x => 0.4728 
K 2 =■ 0.4000 


(such as to produce less than critical damping In each of the two normal 
modes) then the roots of the characteristic equation (16) as solved by 
a computerized polynomial root - finding routine are: 


-6.4019x10 H 

-8.58139xl0“ 4 + j 2.952816xl0 -2 


-1.670404 

verifying the stabiliw/ of the system. 
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V. General Conclusions and Recommendations 

A model is developed for predicting the dynamics of a long, flex- 
ible free— free beam in orbit under the influence of control devices 
which are considered to act at specific points along the beam. Two 
classes of theorems are applied to the system model to establish neces- 
sary and sufficient conditions for controllability depending on whether 
the system possesses non-repeated or repeated eigenvalues. It is 
observed that with a proper selection of the location and number of 
actuators along the beam, a lesser number of actuators than the number 
of modes in the model can control and stabilize the system. 

After establishing the controllability of the system, control gains 
are selected using the following two criteria: (i) decoupling of the 

linearized system equations with appropriate state variable feedback; 
and (ii) applying the linear regulator problem to the modal coordi- 
nates, and thus, selecting gains by solving groups of n n ,f two by two 
matrix Ric^ati equations. 

The decoupling technique avoids modal analysis and is computa- 
tionally simple when the number of actuators is equal to the number 
of modes. However gain selection is possible even when the number of 
actuators is different from the number of modes. The linear regulator 
application described in this report depends on an a priori modal 
analysis and the number of actuators must be equal to the number of 
modes. When the number of actuators is not equal to the number of 
modes the general linear regulator problem can still be applied 
and a 2nx2n matrix Riccati equation has to be solved for a system con- 
taining n modes. 

The Independent modal control concept used earlier for a long 
flexible beam modelled by three discrete masses is reviewed for sta- 
bility when the number of actuators is less than the number of modes. 

For this case, it is seen that even though the system is controllable, 
it is not stable about the zero state vector (gives rise to a simple 

example of control spillover) , It is observed that a proper control ^ 

law not based on modal decoupling ensures stability of all the modes. |j 

l\ 

In the present study control and observation spillover are not \\ 

directly considered and all states are assumed to be available (noise | 

free). Selection of modes for the mathematical model is done on an [ 

arbitrary basis. Only point actuators are modelled. £ 

As an extension of this study, control gam selection using pole 
allocation can be investigated. Model reduction using energy cr 
shape of the structure as a criteria may be studied. Distributed 
actuators can also be modelled and their effectiveness can be compared 
with that of the point actuators. Control and observation spillover 
can be taken into account m designing state estimators and reduced 
order controller designs. 


Appendix A 


Evaluation of Modal Mass (M ): 

n 

The shape function [$ r (x)] of a free-free beam satisfy 

♦J* " x r* r (A’ 1 ) 

, lv * 4 * r 

(where = dx 1 * - and a s ^ m il ar notation denotes other ordered 
derivatives) 

with boundary conditions 


*;'«>) = ^ m (o) 

#”(« - *”•(*) 


(A-2) 


From consideration of Eqs. (A-l) and (A-2), the shape function is 
given by 

(A-3) 

$ r (x) ° cosh X r x + cos A r x - a r (sinh A r x + sin Arx) 
where A ^ is given by the solution of the transcendental equation 

(A-4) 


cos cosh A^fc-1 =» 0 


and 


cosh A r £ - cos \ T i 


sinh A^& — sin A^Jt 

A 6 andV^ tW ° different sha ? e functions 4> r , $ g correspond* n 2 to 
r s 

$ 


(A-5) 


V, 


iJ T - A- 
r r 


(A-6; 

, iv . »* 
s s 

*s 

(A-7) 

* be combined as 


A a IV 

X ,1V 


s v r 

— $ 6 
r y s 

(A-8) 


A-l 






Integrating (A-8) by parts, 

l z 

/ 4> r * s <Ix ■> f Q C'f’a ^r V " ^r *s V]dx 


(X 4 -X i *) *1" - K K" - ^r 4 


r s 


1 1 ? 

r s 


+ K K"h “ 0 for r ^ s (A-9' 

When r = s » n the above integral is defined as nodal mass (M ) per 
unit density per length: n 


M n /P “ / <f>S d X 


(A-10) 


1 , . IV , 

W dx 

n o 


Z Z 

KK” “xl 

n o o 


xti KK" dx 

n o 


Z Z 

!*- i k' c d *] 
" 0 0 


(A-ll) 


4- / <J> * ' 2 dx 


X n n 

n o 


(A-12) 


So 


4M 




(A-13) 


After substitution for i and Its derivatives into Cq. (A-13). 


= p/ dx « pi. 
o 
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Appendix B 


The program described in this appendix solves the equations 
of the form 

X - AX + BU (B-l) 

where 

X is an n dimensional state vector 
A is an nxn matrix 
B is an nxm matrix 
U is a mxl control vector. 

U is obtained using state variable feedback 

i.e. U = GX (B-2) 

using (B-2) , (B-l) can be written as 

X * (A + BG) X 

One can either give (A + BG) as a single matrix to the program or A, 

B, G as separate matrices. The solution is obtained using the state 
transition matrix technique. The plotting is incorporated In the pro- 
gra using the separate computer algorithm REDOK-PLOT (see program 
listing which follows) . 
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6/6/79 15:22:17 

i JOB [READ IN AT 15:20:531 REDOK 

1F0RT ft /8/E/P/S FORT.LS/L ... - 

1LIST1NG 

C. . . . » - FORTRAN STATEMENT 

DIMENSION P!in,101,G!10,101,F!10,101 

COMvON/9EOOY!/PnO,10,iO),NMCC10),ALPHA(lC),RETAnO) ** - • 

01 MENS I ON A(10»!0)«ETGR(t0)»EIGI(10)»CCin#AINV(10,10)»NAME(S) 

Cl MENS ION AK(l0,10)rU(10)»DUMMY(10. 105) ,Z! 1051 , FORCE (5,1 05) 
DIMENSION XflOl,XI!10) 

CALL IN0UTC5/8) 

113 WRITE!*, 11 

1 FORMAT !5X, ' A MATRIX IN THE EOUN. X = A X + RU * ) 

C NiOIMENSION OF 'a' MATRIX ,NP=COLl)MNS OF '5* MATRIX 

r. IF NPCP=1 EQUATION IS 'X=AX' AND NPCP=0 IF X = (A + BG)X' 

C C IF N?c p =l, P and G MATRICES NEED NOT BE GIVEN 

RE AO ( 5, 2 1 N,NP,NPCP 

2 FORMAT f 3 1 2 ) 

WR I TE ( 5 # 3 ) N,NP 

3 FOR w AT f 2X , ' DIMENSION OF A=' , 12, RX, 'COLUMNS OF R = ' ,T2) 

C READING A MATRIX POrtwlSE 

00 a 1=1 ,N 

READ (5, 51 (A(I,J),J=1,N) 

5 FORMA 7( n!r lC.O) 

a WRIT E(«,&) ( A Cl, J) , J = 1 ,N) 

6 FORMAT(2X,!0(P10.a,2X)l 

C READING B MATRIX ROWWISE 

IF (NPCP.GT.O) GOTO 1230 
DO 7 1=1, N 

PE»D(5,5) (SCI, J),J=l,NP) 

7 WP I TE ! 3 » 61 (B(I, J1 , J = 1 ,NP1 - 

C READING G matix OF U=GX ,RQwWTSE 

DO 8 1=1, NP 

REA0(5,51 CG(I,J1,J=1,N) 

8 WRITE!*, 61 !G!I,J1,J=1,N1 

00 9 1=1, N 

OG 9 J= 1 , N 

SUMrO.O 

00 100 m-i,mp 

100 SUM = S'Jm^BCI , M 1*G! M , J1 

9 F!I,„)=S«JN 
DO 11 1=1, N 

11 WRITE!*, 6) (rci,Jl,J=l,N) 

DO 12 I=1,M 

00 12 Jrl ,N 

12 «-! I, Jl=- f I, Jl+AfI,Jl 
WRITE! 3, 1381 

13, F0RMAT!5X , 'MATRIX A=A+PC WHERE U=CX ') 





■ v<* « " 
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DO 120 1=1, N 

20 WRITEC8,61 (AfI,Jl,J=l»N) 

BASIC MATRIX PROGRAMME 

230 CONTINUE . 

NAME OF PROGRAM, ME^*TN 'A' FORMAT 
REA0C5.20M) (NAME Cl), 1 = 1 ,5) 

•001 FOR MAT(SA'j) 

: KEEP A BLANK CARD TO GET ALL Ot-TIONS OF THE PROGRAMME 

: IDET = 0 RRIMTS DETERMINANT' VALUE 

INV=0 PRINTS INVERSE OF A MATRIX 
; NRM=0 PRINTS RES n LVENT MATRIX 

: ICP=0 PRINTS Cl -RCTERISTIC POLYNOMIAL 

: IEIG=0 PRINTS EIGEN VALUE 0 

I ISTM=o PRINTS STATE TRANSITION MATRIX 

: IF ABOVE PARAMETERD ARE NOT ZERO THEN CORRESPONDING VALUES AP£ NOT 

REAO(5,201 3) IDET, INV,NRM, ICP, ICIG, ZSTM 
WR I TE ( 8 * 2008 1 

•008 FORMATrtHloX, * BASIC MATRIX PROGRAM') 
rtRITE(B,200R) (NA M E (I), 1=1,5) 

2000 FORMAT (6X, * PROBLEM TNOENT IFICAT ION : • , 5X , 5 Ail ) 

WRITE (B,2012) 

•012 FORM4T(1HO,4J5(1H*)) 

IF(IOET.NE.O) GOTO 10 
0=0ET( A f N) 

WRITt(8,201 0) 

•010 FORMAT(1HO,SX, * deternimant OF THE MATRIX') 

WRI TE ( 8 » 200 3) D 
•00 3 FoRMAT(1PoE20.7) 
a IF ( IN V.NE.O) GOTO 15 

WRITE(«,201 1 ) 

•Oil FOR M AT ( 1 HQ , 5X , • THE INVERSE OF MATPIX ') 

CALL SIMEQ(A,C,N, AINV,C,IERR) 

IFdERR.E1.01 GOTO 15 
DO 208 1=1, N 

>08 WRITEC8.2003) ( AINV ( I , J) , J=1 ,M) 

.5 CALL CHREK A,N,C,NRM) 

CALL RPOOT(M,C,EIGR,EIGI,+l) 

IF(ICP.NE.O) GOTO 308 
WRITE (8 » 20 l 2 1 • 

WRITEC8, 20051 

>005 FORMA T ( i HO , 5X , ' THE CHAPCTERSTIC POLVNO M I AL-IN ASCENDING POWERP 
1 OF S') 

NN = N+ t 

WR ITE ( A , 20 0 3) (CCn,I-l»NNT 
508 IF(IEIG.NE.O GOTO 35 

«P I TE ( 8 « 20 1 2 J 
wRITE(B,20A6) 

>006 FQRmat ( 1 u 0 , 5X , ' EIGEN VALUES OF A MATRIX ') 

>007 FORMAT ( 9X , 'PEAL PART', 8X, 'IMAGINARY PART') 

>013 ROPMAT(fell) 

WPITE(8,2007) 

00 16 1=1, N 

6 ■vRITE( fl »?0 0 31 EIGPdl ,EIGI (II 

35 TFCI5TM.ME.01 GOTO 25 

CALL ST- , ST(N,A,EIGR,EIGI,IS™1 

B-3 




1 


, +<« 


FORTRAN STATEMENT 


25 

C 

C 

C 

C 

C 

C 

c 

c 

c 

1875 


1724 

1268 


1269 

1270 


30 

40 

80 

20 

1725 

1726 
50 


13 

102 

101 

10 


1271 


CONTINUE 

NsOI«ENSION OF A MATRIX 
MsCCL'JMNS OF 8 MATRIX 
NC=OI*ENSION OF FEEDBACK STATES 

NwRITEsO IF PRINTING IS NEEDED * * 

NPLOTsO IF PLOTTING IS NEEDED 
NSEI.CTsO ALWAYS 
TsINITIAL time 

t«ax=final time ~ * - — . . _ 

Xuin-mjnimum or SUM OF states 

read C5, 1875) N,m,nc,NWRITE,NPLOT,NSFLCT 

F0RMATC2I2, 14,311)- - - - 

READ{S,20HXl(I),I = t,N) 

READC5,20) T, Tmax,h, xmjn 

IF(NSELCT.GT.O) GOTO 1724 - • - 

READ (5, 20) (ALPHA(I),BETA(I) ,I=l,N) 

DO 1268 1=1, M 

READ(5,20r CAKfr,J),Jsi;M) 

WRITE(R,40) CXI(I),Isi,M) 
wRITE(8»40) T ,TMAX,H,Xmjn 
WRITE (8,40) ( ALPHA (I),8ETA (I), 1 = 1 ,N> 
wRITECfl, 1 270) 

00 1269 J = 1 , N 

WRITE (8, 40) (AK(I, J) ,I=l,v) 

FORMAT f 2X » 'K MATRIX WRITTEN COLUMN WISE') 
IF(NSELCT.GT.O) GOTO 1 725 
RE AO C 5, 80 ) CNVC(I) ,1 = ) ,N) 

00 30 1=1, N 
00 30 J=1 ,M 

RF AO (5, 20 ) (P(I,J,K),K=1,N) 

F0RmaT(2X, 1P6E20.7) 

F 0 ° v 4 T (301 1 ) 

F '°ma T ( 0F 1 0 . 0) 

DO 1726 Nmk=1,N 
DO 1726 Nml=1,N 

WRITE (8, 40) (P(NMK,N m L,Nmm) ,N v > >1=1 ,N) 

N m N r 1 

00 10 J= 1 , N 

X(J)=0.0 - - . 

DO 101 1=1, N 
A0P=4LPHA(I)*T 

IF(A5P,GT. 100.0. OR. ABP.L-T. -100.0) GOTO 101 
DO 102 K = 1 , N 
IF(NMCCI) .EO.O) GOTO 13 

XCJ)=XCJ)+P(l,j,K)*CXPCALPHA(I)*T)*C0S(9£TACf)<.T)*XICK) 
GO TO 1 02 

XCJ)=X(J)+P(I,J,K)*EXPCALPHA(I)*T)*SIN(RETA(I)*T)*XI(K) 
continue ■ * — 

C 0 N T 1 N U E 
CONTINUE 
00 1271 1=1, m 

U(I)=0.0 

DO 1271 J=1,N 

U(I)=X(J)*AK(I,J)+U(I) 

00 1371 K = 1 , m 


•n 

u 
\ 


5 

V 


% 


i 




M 

5 
i 
t 
i 

t 


A 

*1 






< * 


C .+< FORTRAN STATEMENT — — N 

* i 

1371 FORCE(K,NMN)=U(K; 

IF(N.nRITE.GT.O) GOTO 1000 
WRITEC8»t03) T, (X(J),J=1»N} , OKI )> 1 = 1 » M ) 

103 FCRMATC2X, 1P8E16.5) 
lflOO CONTINUE 

DO 1272 1=1, N ' “ 

1272 DUMMY(I,N^N)=X(I) 

SUM=0.0 

DO 70 «w=l,N * ' *" * 

70 SUM=SUN+X 
T=T + H 
NM,M = NMN*l 

Ir(SU v '.LT.X M IN.OR,T.GT . T M AX) GOTO 60 
GOTO 50 

60 CONTINUE 

CALL FOPENd , "DPOtREOOK") 

WRITE 9INAPY(nN,NC»M 

.iRITE BINARY ( 1 ) C ( DU M MY ( I , J ) , J= l , NC ) , I = 1 , N ) 

WRITE OINARYCn ((FORCECT,J),J=l, N C),I=l, M ) 

CALL FCLOSE(l) 

GOTO 11 3 

1020 STOP _ . 

END . . 
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C. . ..+<-** FORTRAN STATEMENT — _-> 

SUBROUTINE CHREO(A,N,C,NPM) 

C THIS SUBROUTINE FINOS THE COEFFICIENTS 0* THE CHARACTERISTIC POLY 
C NO M I AL USING THE LEVERRIER ALGORITHM 

COMMON ZED (10, W*, 10) 

DIMENSION A!10,10),C!11),ATEMP!10,10),PROD!10,10) • 

•1000 FORMAT! 1H0,5X, 'THE MATRIX CCEFFICENTS OF THE NUMERATOROF THE PESO' 
1LVENT MATRIX ') 

1001 FORMAT! 1H0,5X, 'THE MATRIX COEFFICENTS OF S', II/) 

10*2 FORMA f f 1P6E20.7) - * 

1003 FORMAT! lH0,a5(lH*) ) 

C REPLACING THR DATA CARD DATA A TEMP/ 1 00*0 . 0/ 

DO 1315 in, 1C - ----- -- - - * 

DO 1315 J = 1 - 10 
1315 ATEmpti, o 0 

CALL CHM . , N,C) - - - 

DO 65 I = « - 

65 ATEMPTI,. : o 

70 00 80 1=1 , U 
00 80 J=1,N 

80 ZED(N,I,J)=A-TE m P(I,J) 

IF!NRM.fjE.O) GOTO 71 
WRITE!*, 100 3) 

WRITE!*, 10P0) 

M = N >1 

WRITE!*, 1001 ) M 
DO 35 1=1, N 

35 WRITE!8, ICO?) ! ATEMPIII, , J) Y9J=1 ,N) 

SYNTAX IN ERROR, PU v CTUATION MISSING, OR IDENTIFIER OF WRONG VARIETY 

71 DO 40 1=1, N 

DO a 0 J = 1,N 

ao ATE V P!I,J)=A!I,J) 

DO 10 1=1 , N 
NNN=N-I 

IFd.EO.n GOTO 55 
IFINRM.NE.O) GOTO 60 

WRITEI8, 1001 ) NNN - - - - - - 

DO U5 J=1,N 

«5 WRITEI8, 100?) ! ATEMP!J,K) ,K=1 ,N) 

60 NP=NNN + 1 • - 

DO 90 11=1, N 

00 90 J = 1 , N 

90 ZEDINP, II, J)=ATEMR!II,J) 

DO 15 J = 1 , N 
DO 15 * = 1,n 
PROD! J,K)=0.0 
DO 15 L=1 , N 

15 PROD! J,K)=PP00! J,K) + !A! J,L)*ATEmp( L ,k) ) 

DO 1 3 J= 1 » N 
DO l 3 K=1,N 

13 ATEMP! J,K)=PPOD! J,K) 

55 00 10 J=1,h 

10 A TEMP! J, J)=ATEMP!J, J1+C!N-Itn 






> 


* 



SUBROUTINE CHMRQA£A,N,C) 

DIMENSION JCll),C(ll),8(10,10),A(lO,10)»Df 

NN=N+1 

00 20 1=1, nn 

0 C(I)=0.0 

C (NN) = 1.0 

DO la M=t,N 
K = 0 

J ( 1 ) -1 

GOTO 2 

-JCUsJCLUI ' 

• .IFCL-M) 3,5,50 

3 MM=M-1 

DO U I=L, M M 

II=I+i 

1 J(II)=J(I)«-1 

; do 10 i=i, M 

00 10 KK=1,M 
NR=J C 1 5 

- NC=J£KK) * 

.0 8CI,KK)=A(NP,NC) 

K=K + 1 

D (K 1 =D£T ( 0 , M ) 

DO o 1=1, M 
tsM-I*| 

IF(JCL)-(N-M+L) ) 1,6,50 - 
, CONTIN UE 

1 

do i« i=i, * • ‘ 

1« C( M l)=C( M l)+D(I)*(-l.O)**M 

RETURN 

50 WRITE (8, 2 000) 

>000 FOR M ATC 1H0 , 5X , * ERROR IN CHREQA ) 

RETURN 

• - END - — 


300 ) 




FORTRAN statement 


T«s”S''!c?I^ # SOTPRO G BAM FINOS THE 06TER-IN.NT OF <_WI> 
USING oiagonalisation procedure 
DIMENSION Alio, 10 ), 8(10,10) 

IREVsO , • T 

DO 1 1=1, KC 

DO 1 J=l,KC 

B(I,J)=A(I,J) 

• DO 20 1 = 1, KC T ' ' : * 

K=I 

IF(8(K,I)) 10,11,10 .. - 

• K-=K + 1 

IF(K-KC) 9,9,51 

IFCI-K) 12,1^,51 

DO 1 3 m= 1 , KC 

TEMP=8CI,M) " 

B(I,M)=0(K,M) 

B(K',M)=TE v P • 

IREV=IREV+1 

11=1+1 . ~ .. 

IF(II.GT.KC) GOTO 20 ~~ * 

DO 17 M=I I , KC 

IF (B(M f I) ) 10,17,19 _ 

TEMP=0(M,I)/B(I,I) 

DO 16 N=I,KC 

8(M,N)=0(M,N)-B(I/N)*TE V P _ . . 

CONTINUE 

CONTINUE 

DET=1 . - • 

DO 2 1=1, KC 
DET=DET*B(I,I) 

DET=(-1 )**IREV*DET . . . 

return 

DET=0.0 

return _ „ . 

• END ~ " * ~~ ' ~ 
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SUBROUTINE PROOTIN, A,U,V,IR) nnnTO 

THIS SUBROUTINE USES A ^OOIFIEO BARSTO^ METHOD TO FIMO THE OOTS 

OF A POLYNOMIAL. 

DIMENSION A(20),U(20)»V(20)»H(2l), 9(211,0(21). 

IREVsIR ** 

NC=N+1 

00 1 1=1, NC 
H(I)=A(I) 

PrO.O “ * 

0 = 0 . 

R = 0. 

IF(H(m a, 2, a 
NC=NC-1 
V (NC ) =0 , 

U(NC)=0. 

00 1002 1=1, NC 
H(i)=H(i+n 
GOTO 3 

IF(NC-l) 5,100,5 
IF (NC-2) 7,6,7 
R=-H(l)/H(2) 

GOTO 50 

IFCNC-3) R,e,9 

P=H (2 ) /H ( 3) 

0=H( 1 ) /H(3) 

GOTO 70 

IF(A8SCH(NC-n/H(NC))-ABS(H(2)/H(l))) 10,19,19 

If?EV = -IREV 

M=NC/2 

DO 11 1=1, w 

NL-NC-I+l 

F=H(NL) 

H(.NL)=HCI) 

H ( I ) =F 

IF (Q ) 13,12,13 
P = 0. 

GOTO 15 
P=P/Q 
0=1. /O 

IF(R) 16,19,16 
R=l./R 

E=5.E-10 

BfNC)=H(MC) 

C(NC)=H(NC1 

3(NC+1)=0. 

C (NC+ l ) =0 . 

NP=NC-l 

no r 9 j=i,iooo 

00 21 11=1, NP 
I='IC-1 1 

B(I)=H(I)+R* R (H 1 ) 

CCI)=9(1UR*C(IM) 

IF ( ABS C R ( 1) / H (l ) )-E) 50,50,2^ 

I F C C ( 2 ) ) 23,22,23 
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• + < __ FORTRAN STATE w cNT 


>2 

>3 

30 


37 

* 3t 

32 

33 
3« 


35 


36 

OR 


50 

51 

52 

53 
5« 

70 

71 

72 

73 
7a 


75 

80 


81 

82 


76 

77 


R=R + 1 

GOTO 30 __ 

RsR-9 ( 1 ) /C ( 2 )“ ‘ 

00 37 IlslfNP 
I=NC-Il 

8 ( I) =H ( I) -P *B ( I + ( 1+2) 

C(I)=9(I)- P *C(I+1 )-0*CCI+2) 
IF:h( 2)) 32,31,32 
IFt A8S(8(2) /H( 1 ) )-E) 33,33, 
IF(A8S(0C2)/H(2))-E) 33,33,3a 
IF(A8S(9(1)/H(1))-E) 70,70#3a 
C8AR=C(2)-9(2) 

0=CC3)**2-CBAR*C(a) 

IF (0) 36,35, 36 

P=P- 2 

Q=G*(QM .) 

GOTO a 9 

P=P+(8(2)*C(3)-B(l)*CCa))/0 
O=Q+(-8C2)*CBAR+0Cl)*C(3))/O 
CONTINUE 
E = EMO. 

GOTO 20 
NC=NC-1 
V(NC)=0. 

IFCIREV) 51 ,52,52 
U(NC)=1 ,/R 
GOTO 53 
U(NC)=R 
oo 5a 1 = 1 , nc 
H d)=9(i + n 
GOTO a 
NC=NC-2 

IFCIREV) 71,72,72 

0P=1 ./O 

PP = ° / ( Q*2 .0 ) 

GOTO 73 
QP=Q 


FP=P/2. _ 

F=(PP) **2-QP— 

IFCF) 7a, 75, 75 

U (NC+ 1 ) =»PP __ 

U(NC)=-P° ' 

V(NC+1)=SQRT(-F) 

V(NC)s-V(NC+l) _ . 

GOTO 76 

IFCPP) «1,80,P1 

U(NC+1)=-S5RT ( F) 

U(NC+1 )=-(PP/A6S (PP) )* ( ABS(PP) + SQRT (F) ) 

CONTINUE 

v(Nc+n=o. 

U(NC)=QP/U(NC+1 ) 


V (NC) =0 . 

00 77 1=1, NC 
H(I )=8CI+2) 
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GOTO a 
00 RETURN 
END 
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FORTRAN STATEMENT 


SUBROUTINE STMST(N#A,EIGR»EI6T,IKN0W) MATRIX USING 

THIS SUBROUTINE DETERMINES THE STATE TRANSITON MATRIX USING _ 

SYLVESTER'S EXPANSION THEORE M 

COMMON CHI ( 10, 10,10) oi:TI , M(n 

COMM0N/REOOY 1/P C 10^ 0 , 1 0) ,NMC C 1 0) , ALPHA ( 10) , BET ,( . 

no) - 

11PE13.6,')T»)- 

III MATRIX COEFFICIENT OFEXF ( ■ * 'Jl 

)0» 'FOR-tuiHO^EX, -THE MATRIX COEFFICIENT OF EXPt ‘ , 1PEI3.6, •>! ■ ) 

JOS FOR M AT(lHO,ii5(lH*)) 

IFUKNOW.NE ,0) GOTO 800 
WRITE(8,100S) 

JO 00 10 K = t,N 

CEIG (K ) sC^PLX (EIGR(K) ,EIGI (K) ) 

00 10 L=1,N 

D CA(K,L)=C M PLX(A(K,L),0.0) 

1 = 1 

IF(IKNOW.NE.O) GOTO 700 
WRITEC8, 1000) 

00 00 15 «=1,N 

5 OENOM(K)=CeiG(I)-CEIG(K) 

00 500"J=l,N 
IF(J-I) 1 00,500,200 
00 IF(J-t) 1*0,110,150 

00 IF(I-l) 300 , 300 , aOO 

00 IF(J-I-t) 110,110,150 

00 IFCJ-I-1) 110,150,150 

10 DO 5 K = 1,N _ .. . 

DO 5 L=t , N 

CA1 (K,L)=CA(K,L) 

DO 20 K = 1 , N -- - ~ - - - - — 

CA1 (K,K)=CA(K-,X)-CEIG(J) 

DO 20 L=1,N 

0 CA1 (K,L)=CA1 (K,L)/0EN0M(J) 

GOTO 500 

50 DO 00 K=1,N 

DO O0 L = 1,N _ 

0 rA?(K,L)=CA(K,L) 

DO 25 K=1,N 

CA2(K,K)=CA (K,K)-CEIG(J) 

DO 25 L=!,N 

5 CA2(K,L)=CA2(K,L)/DENOM(J) 

DO 30 *=1,N 
DO 30 L = 1 » N 
TCA(K,L)=(0. 0,0.0) 

00 30 M =1 9 M 

,0 TCA(<,L)=TCu( l ',L)+CAl(K,M)*CA2(°,L) 


1 
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. .+< 

DO 35 K = 1,N 

00 35 L=1 ,N 

35 CA1 (K,L)=TC4(K,L) 

500 CONTINUE 

IF(AImag(CEIG(I))) <15, SO, 05 

05 IM=i 

1 = 1 + 1 

ALPHA(IMN)=EIGR(I) 

* ~BETA(IMN)rEIGI (I) * 

NMC(IMN)=1 

IF(IKNOW.NE.O) GOTO 80t 
WRITE(3,1001) EIGR(I)fEIGICI) 
801 00 65 K = 1,N 

DO 65 L=1.N 

65 SPS(K,U=REAL(CA1(K,U)*2.0 
00 66 K=t ,N 

00 66 1=1 , N 
CHI(IN,k,L)=SPS(K,L) 

66 CONTINUE 

00 1100 J=1,M 

00 1100 * = 1,N 
1100 Pfl^N, J,K)=S P S(J,K) 

IM^rjMfj+1 

IF(IKNOW.NE.O) GOTO A02 
00 80 K=1 ,N 

80 WRITEI8, 1002) ( SPS (K ,L) ,L= 1 , N) 

ALPHA(IVN)=EIGR(I) 
BETA(IWN)=EIGIfI) 

NMCIIMMjsO 

WRITEC8, 1003) E IGR C I) , EIGI ( I ) 
B02 DO 55 K = 1,N 

00 55 1 = 1, N 

55 ~ S P S(K,L)=AI M AG(CA1 (K,L))*2.0 

00 56 Ks 1 , N 
00 56 1 = 1, N 

CHICI,K,L)=SPS(K,L) - 

36 CONTINUE 

OC 1110 J = 1,N 

00 1110 K = 1 , M 

1110 PCIMfj, J,K)=SPS(J,K) 

JMNsIMN+1 

IF(IKNOW.NE.O) GOTO 600 
DO 85 K=1,N 

15 W p ITE(8,1002) (SPS(K,L),L=t,N) 

GOTO 600 
.0 CONTINUE 

ALPHA(IMN )=EIC ,R ( n 
BETA ( I M N) sO . 0 
N«C(INN)sl 

IF(IKNOW.NE.O) GOTO 300 
WRITE (8, 1004) EIGR(I) 

00 00 60 K = l,to 

00 60 L = 1 t N 

0 S p S(K,L)=REALCCAl (K,L) ) 

DO 61 K = 1 , N 
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C. . . FORTRAN STATEMENT 

DO 61 L=1 r M 
CMI(I»K f L)aSPS(K,L) 

61 CONTINUE ' ' - 

DO 1120 J=lfN 
00 1120 K=1,N -* 

1120 P(I.MN r J,K)=SPS(J,K) ’ 

IMNsIMNM 

IF(IKNOW.NE.O) GOTO 600 

— DO 75 Ksl,M 

75 WRITECS, 1002) (SPS (K ,L) ,L=1 » N) 

600 TFCI.GE.N) RETURN 

-1 = 1 + 1 *’ 

GOTO 700 
END 
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C....+< FOPTRM STATEMENT 

SUBROUTINE SI w EQ(A,XDOT,KC, AIN V,X, TERR) 

C THIS SUBROUTINE t INDS THE INVERSE Or THE MATRIX A USING 

C OIAGONALI2ATION procedures 

01 MENS I ON A ( 10, I 0) ,BC10, 10) ,XOOT Cll) ,X C 11 ) , AINV(JQ,10) 
Nr 1 -** 

IERR=1 
DO 1 1=1, KC 

00 l J=l,KC __ 

AINV(I,J)=Ov * 

1 8CI,J)=A(I,J) • 

DO 2 1=1, KC 

AIMV(I,I)rl. 

2 X(I)rXnOTCI) 

00 3 1 = 1, KC ’ 

-• " - COMprO. 

K = I 

IF(A5S(8(K,I))-ABS(C0 N| P)) 5,5,0 

comps0(k,i) 

NrK 
K = K + 1 

IF(K-KC) 6,6,/ 

IF (B(N, I ) ) 8, SI, 8 
IF(N-I) 51,12,9 
DO 10 m=i,XC 
TE M P=B C I , w ) 

6CI,M)rB(N,V) 

B(N,M)=TE M P 
TEMP=AINV( I, M ) 

AINV(I,M)=AINV(U,M) 

JO A I N V ( N , M ) = T E MR 

TEMPrX(I) 

X(I)=XfN) 

X (N) =TE M P 

12 X(I)=X(I)/8(I,I) 

TEMP=B(I,I) 

DO 13 M =1,KC 

A I N V C I , w ) = A I N V f I , M ) /TEMP 

13 3CI,M)=BCI,mi/TEMP 

DO 16 J = l ,>-'C • " 

IF(J-I) 10,16,10 
10 IF(8(J,D) 15,16,15 

15 X(J)-X(J)-3(J, * ) *X ( I ) 

TE^PrB ( J , I ) 

DO 17 '» = 1 , K C 

AiMV(J,M=AINV(J,N)-TE^P*MNV(I,N) - 

17 8CJ,N)=9CJ,N)-TE^ *CCI,N) 

16 CONTINUE 

3 CONTINUE 

RETURN 

51 -r:TE(»,52) 

52 FOR*'A TC6X, 'THE MATRIX IS SINGULAP ') 

IFRPrO 
RETU°'l 
FNO 
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SLHUUL CJF ENGINEERING 


*******.»****»•*•******************«***********•**#****»********************■*»• 

2/26/79 1 2: 12;22 ' 

! J08~ CRF 40 IM i’T 1 2 * Q t 33] "R£D0 k~PLOT PROGPAM 

IFOR T /A /9/E/P/S FORT.LS/L . . 

•LISTING ** 




FORTRAN STATEMENT 


IN THE DATA CARO /READ STATEMENT NA DECIDES BEGINING OF X AND XOOT 
CURVES , Nr 0 SKIPS ALL X-CURVES ma DECIDES BEGINING OF FORCE CURVES 
SKIPS ALL FOOCE CURVES. IF N.NE.O AND m.NE./ DECIDES HOW many FO 
DEFLECTION CURVES THEY NEED 
DI M ENSION dummy < 10, t OS), FORCE C5, 105) , 2 (105) 

CALL I NOU T ( 2 , 8) 

CALL FOPEN ( 1 V"DP0 : PEDQK " ) 

READ BINARY Cl) NX,NC,MX 

READ BINARY ( 1 ) ((DIJWYCI, J),J=l,NC),I=t,NX) 

READ 81 NAR /Cl) ((FORCECI»J)f J=1»NC),I=1, M X) 

CALL FCL05ECU 

RE AO ( 2 » 26 7S ) NA,N,m A ,m 

FORMA T ( il I 2) 

IF(M.EO.O) GOTO 2676 
DO 127a 1=1, N 
DO 1 275 ~J = 1 , NC 
ZCJ)=0'JMMY(I,J] 

ZMAX=Z(J) 

ZMINsZd) * 

00 1276 IJ=lfNC 

IFCZMAX.LT.Z(IJ+1)) GOTO 1277 
GOTO 129 s 
ZNAXrZ(IJ+l) 

IF(Z M IN.LT.Z(TJ+11) GOTO 1276 
Z'*I\ = Z ( I J + 1 ) 

CONTINUE 

IF ( A8S (ZMAQ , GT. ABSCZMINn GOTO 6376 

ZMAX=A9S(ZMIM) 

GOTO 6379 

Z m X M r - A G S_( Z M A X 1 

CONTINUE 

WRITE («,15ni) Z W AX , Z v I N 

CALL PGR ID ( n_,J? ) 

CALL P L 0 X ( Z M I n", Z , Z M A x , N C ) 

CALL °LOGO (0,0,6.01 

.CONTINUE 

CONTINUE' 

IF C V .EO.O) GOTO iaio 
on 1 376 1=1 , u 
DO 1377 J = ! >C 
ZCJ)=FORCE(I,J) 

ZMAx=zm_ 

Z m I n = z c t ) ' 



--- FORTRAN STATEMENT 

0° H73 T.TK = i-,vr 



— 7^AX= z c I_JK <- n 

iF(ZMir) tL r7?TnK7IT, ~:c , 37 q L 

ZMINsZf IjK+l) ' 578 

CONT TN IIP 

goto 635o — 

— G OTO fe ^qt 

ZM I fJ =-a^s(ZMATj— ■ — 

CONTINUE 

150n Z M AX . Z M I H 

f orm* r (~X, — 

C^LL PGRIQCO,?) '" 2 °* 7 ' 3X, »**z«I.V=* # E? 

-f ff r-P- 9 * c -? M I^>Z > ZM flx . Nr ) 

CALL PLOGQCit70,6.0) — 

continue 

_CONTINijE 

CALL EXIT 


Appendix C 

The program described in this appendix solves the non-linear 
equations of motion of the beam in orbit incorporating the control 
laws obtained using the linearized model. 

It can plot deflection of the beam at various instants of 
time with control and the time history of actuator forces required. 

Data cards to be given are explained in the program by means 
of comment cards. A listing of this program follows. 







HOWARD UNIVERSITY -- SCHOOL OF ENGINEERING — 
***********»**************************************•***********************' 
6/7" t S : 3 3 S 8 ~ 

'08 tREAO IN AT 15:32:16) NON-LINEAR 

0RT/A/8/E/P/S FOPT.LS/L - — - * - - 

ISTING 

. . .♦< — --FOPTRAN STATEMENT > — 

external FCT,OUTP 

COMHON/REOO/AOMEGACIO) *OAMPC10),WC # NMODES»LNEAP,NACT,MSELCT, 1ST ate 
COMMON/RED/Z(10,20) 

COMHON/REDOO/FREQCIO) ,Q(?''),NP0NTS,ALC20),ISELCT,NPL0T,NWRITE 

COMMON/RAJ/QMAX 1 ,QMIN1 - • 

COM^ON/REAM/PTTCHt 100) ,NPK,QS«ALL 
COMMqn/REJA/AK (5# 10)#F(5#100)#FF(100) 

common/ki^/imult ' • - - • 

DIMENSION YC?0),OERY(?0),AtlX(e,20),A(a) f BCU),CCa),PRMT(5) 

DIMENSION SIZEC10) 

CALL IN0UT(5,«) 

<11 CONTINUE 

CALL PLOGOCO. 0,6.0) 

Q M AX 1=0.0 
NPK = 0 
QMINlsO.O 

WC=0R5ITAL FROUENCY,TOL=TOLRERANCE FOP RUNGE-*UTTA SUBROUTINE# 
QSMALL=SM4LLEST VALUE OF DEFLECTION ALONG THE BEA" TO STOP PLOTTIN 
SIZECD-MAXIVIJM VALUES OF STATES FOR RIJNGE-KUTTA ROUTINE 
AOMEGA(I) FRQUENCY VALUES 

PRMTC l )=INm al TIME, PRMT(2)=FINAL T Im£ , PRmt ( 3 ) =1 NCRE^ENT 
Z(M#N)=FEEOBACK gain matrix 
• FREQ1I) FREQUENCY VALUES" • 

AL(I)_ "POSITIONS ALONG THE BEAM at WHICH DEFLECTION IS CALUCULATE 

readisTrooo) wc#tol,osmall 

NSELCTsO IF DEFLECTION PLOTS ARE NEEDED AT ONE PLACE OTHEPwlSE ONE 

N MOD ES= NUMBER OF MODES CONSIDERED INCLUDING PITCH 

NPONTS = NiJMRER OF POINTS-ALOMG THE- BEAM 

I SELCT = T I M E INTERVAL SELCTION OF PLOTING 

NPL0T=1 IF PLOTTING IS NEEDED OTHERWISE ZERO 

N.nRITEsI IF WRITING IS NEEOED OTHERWISE ZERO 

NLNEARsl IF EQUATIONS ARE NON LINEAR OTHERWISE 7ER0 

QSMALL=QUANTTTY DEFINING the SMALLEST maximum DEFLECTION ONE WANTS 
TO PLOT (DEFINING PRACTICAL ZERO) 

NACT = N'JMRFR OF ACTUATORS 

NACT=Q THEN DATA CARDS FOR AK(I,J) NEED NOT RE SUPPLIED 
CONTROL FOPCES CAN MOT 8E PLOTTED IF NACT=0 
ISTATEsNO. OF. STATES CONSIOERED FOR FEEDBACK 

RE ADC 5,00 la) NmqDES . NPONTS , I SELCT , NP1.0T , N.tO I TE , ML^E A», NACT , NSFLC T 
1 rTSTATF, ImijlT 
I«k= 2 frNMPDES 
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enoTDAM ptatp.mpmT 


> 



9000 FORMAT ( 3F 1 0 . 0 ) 

9003 F0R M AT(2X,I2»2(E13.6,2X)/) 

9009 F0RMAT(?X,6(E13.6,2X)/) 

9005 F0R M AT(50X, 'SIZE VALUES') 

9006 FOR^ATCSOX, 'FR&UENCIES') 

9007- FOR M ATC20X,-'PAPAMETERS:INITIAL TI^E, FINAL T IME , INTERVAL ' ) 
9009 FOR M ATC50X, 'IMITIAL VALUES') 

9011 FORMATC20X, 'Z-BOTTO^ PART 0 e A-MATRIX') 

9123 * FORMAT (SOX, 'NUMBER OF BISECT IONS *, 1 2 ) 

9001 FORMAT { 8F 1 0 . 0 ) 

REAO(5,900t) (SIZECI),I=1,IKK) 

READ (5, 900 1 ) (AO y EGAfI),I = l,NMQOES) 

REA0(5,9001) CPRMT(I),I=1,3) 

READ(5,900t) (YU), 1 = 1, IKK) 

CO 108 M3l v N»ODES “ 

108 PEAO(5,9001) (Z(W,N),N=1,IKK) 

9019 FORMATC2I2, 19,612, la) 

- • PEA0C5, 9001) - (FREQ(I), 1=1, M^OOES) 

READ(5,9001) (AL(I),I=l,NPONTS) 

WRITE(9,903t) 

9031 FORMATflOX, 'VALUES OF K MATRIX') 

IF(NACT.EO.O) GOTO 9037 
DO 9030 I=1,m&CT 

RE AO (5, 9001) (AK(I,J),J=1,ISTATE) 

9030 WRIT£(9,9009) ( A K ( I , J) , J = 1 , I ST A TE ) 

9037 CONTINUE 

WRITE! 9, 9003) NMODES,WC,TOL 
WRITE(fl,9O05) 

WRITE (?» 9009 ) (SIZE (I), 1 = 1, IKK) 

WR ITE(9, Q 006) 

wRITE( 9,900<4) (AO*^EGACI),I = l ,N m OOES) 

WRITE(8,9007) 

wRITE(9,°009) (PP M T (I ) , 1=1 ,3) 

WRITE (9, 9Q09) 

WRITE (0,9009) (YU), 1 = 1, IKK) 

WRITE(8,901 1 ) 

00 9013 M = 1 , NMQDES 
WRITE(8,9009) (ZC M ,N),N=1, IKK) 

9013 CONTINUE 

IF(NSELCT.EO.t) GOTO 9092 
CALL PGRI 0 (0,2) 

9092 CONTINUE 

CALL PKSCL(IKK,SIZE,OERY,TOL,PR m T) 

CALL RKGS(PR'-‘T, Y,OERY, IKK , I “LF , FC T , OUTP , AUX ) 

PVAX=PITCW( 1 ) 

PVTNrPI TCH( 1 ) 

DO 0015 IL = 2 » NPK 
IF (PMAV.LT.PITCW(IL) ) GOTO 9016 
GOTO 9017 

9016 P v AXrPITCH( IL) 

9017 IF(Pmin # lT.PITCw(IL) ) GOTO 9015 
PMINrPI TCM( IL) 

9015 CON riN'JF 

IF(A = S( 3M AV) ,C-T. A3S( S,M IN) ) GOTO 901* 

PMAX=A0SfP M IN) 
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GOTO 9019 

9018 PMIN--p\»ax 

9019 CONTINUE ..... _ 

WRITE(8,9020) P«M,PMIN 

9020 FORMAT ft OX, »**PMAX=»,E13.6,5X,'**PMIN=»,E13.6) 

IF (NSEICT .EG.! ) G0T0-90tt3'-' 

CALL PLOGOCO.O, a. 75) 

)0«3 CALL PLOGOCO.O, 1.25) 

CALL PGR 10 ( 0,2)' --- - -- _____ 

CALL PLOX(PMIM, PITCH, PMAX,NPK) 

WRITE(8,9123) IHLF 

FFMAXlsO.O . .... 

FF M I N 1=0.0 

IF (NACT.E9.0) GOTO 9036 

DO 9021 1 = 1 , N ACT -- _. . 

DO 9022 J=1,NPK 
'022 FF (J)=F(I,J) 

FFMAXrFF(t) • - - - - - . _ 

FFMIN=FF(1) 

DO 9023 11=1, NPK 
IF{FFMAX,LT. FF(II+l)) GOTO 902a 
GOTO 9025 

029 FFMAX=FF(H + 1 ) 

025 IF(FFMIN.LT.FF(II+!) ) GOTO 9023 
FFMINsFF ( 1 1 + 1 ) 

023 CONTINUE 

IF( ABS(FFMay) .GT. ARSCFFMIN) ) GOTO 9026 
FFMAX=A8S(FFMIM) 

GOTO 9027 

026 FFMINr-Af3S\.FF M AX) 

027 CONTINUE 

WRITE ( 8,9035) FFMAX,FFMIn 

035 FOR^ATCEX, ' **FFMftX= * ,E 1 3 . 6, 5X , '**FFMIN= » , El 3.6) - 

IF (NSCLCT .EQ.l ) GOTO 90aa 
CALL PLOGOCO.O, a. 75) 

Oaa CALL PLOGOCO.O, 1 .25) 

CALL PGRI 0 C 0 , 2 ) 

CALL PLOX CrFMIN , F ~ , FFMAX , NPK ) 

021 CONTINUE - --- . _ 

WRITEC8,R033) 

D33 FORMAT C2X, 'VALUES OF FORCES STARTING AT TMIE = 0.0') 

DO 90 32 J = 1,MPk 

)32 WRITEC8,900a) ( F C I , J) , 1=1 , NACT ) 

)36 CONTINUE 
GOTO 9001 
TOO CONTINUE 
STOP 
END 

ISTATE=NO. OF. STATES CONSIDERED FOR FEEDBACK 

REA0C5.90 1 U) Nv«oDES,NPONTS, ISELCT,NPLOT,N.vRITE,NLMEAR,MACT.NSPl FT 
l.ISTATE, I M, JL T L 

IKK=,?*NV00ES 
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. ■ FORTRAN STATEMENT -- — -> 

SUBROUTINE FCT (X , Y , DERY ) 

COMMON /REDD /AOMEJjjA ( 1 0 ) , DAMP ( j o ) , WC , N^OOES »LNE AR , MAC T , MSELCT 
CnMM 0 N/REO/Zno,20) 

COMMON/REJA/A* C5, 10),FCS,100),FF(100) 

COMMON/KIM/I^ult 

DIMENSION Y(20}*DERY(20)fAUX(8»?0)»A(U) f B(0) # C(C),PRMT(5) 

IF(NLNEAR.EQ.O) GOTO 10 

Z(lrl)s-1.?*SIN(2.*Y(!)) 

00 10 I =2 t M M ODE$ 

Z(I»I)=-CCA0^EGA(in**2/(WC**2)-(3.*CCSINCY(n ))** 2 )-l ,)-((Y(NMOOF 
is+n/wc-n** 2 n 

10 CONTINUE 
GOTO 15 

a Z C 1 # 1 ) s»3. 0*Y Cl) 

00 16 Is2,MV0t)ES 
Z(I» I )=-( AOMEGA(I) )**2/(WC**2) 

CONTINUE 

DO 11 1=1 ,N v ODES 

1 DERY(I)=Y(I+NMODES) 

PITCH EQUATION 

NKC = N'«00ES+1 

0ERY(NKC)=Z(1 , n+Z(l,NKC)*Y(NKC) 

00 1? I=2,NM0PES 
Ijsl+NVOOES 
IK=2*N M 00ES 
DEPY(I J)=0.0 
00 13 Jrl,IK 

OEPYd J)=DERY(IJ)+ZCI,J)*Y(J) 

3 CONTINUE 

2 CONTINUE 
RETURN 

END 
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SUBROUTINE OUTPCX, Y,DERY,IHLF,NniM,PRMT) 

LOGICAL Rk'NXT 

CO MM ON/REOD/AO w EGA (10)»DAMP(10),WC# NICHES rLNE AP , NACT » NSELCT , 1ST ATE 

COMmom/P£ 0 do/FPEQ(10),OC20),NPONTS, *L(?0) , ISELCT,NPLOT#NWRITE 

CQMMON/RAJ/QMAYl^rjMINl 

COMMON/BE Am/P ITCH (1 00), NPK,QSMALL 

COMMon/r E j A / ak (5, 1 0 ),F (5,1 00 ),FF Cl 00) 

COMMON /k I m/ I m(JLT 

OIMENSION Y(20D ,OERY(?0»AIJXC8,20)-, A(a)-,B(U),C(«)rPRMT(5) 

IF(.NOT.RKNXTCIHLF)) GOTO 2 
IF(NWRITE.EO.O) GOTO 13 

WRITEC9,R01) X, (Y(I) ,I=l,NOIM) ’ ‘ 

01 F0RMAT(2X,F10.3,6CE13.6,2X)) 

3 CONTINUE 

I JK=2*NMQDES - - 

ISSsISTATE/2 
IF(NPLOT.EQ.O) GOTO 12 

IMULTsIMULT+1 - * * 

IFdMULT.EO.n GOTO 11 
IX=IMULT/ISELCT 
IMsIX*ISELCT 
IF(IM.EG.IMULT) GOTO 11 
GOTO 1? 

11 CONTINUE 

NPKsNPK+l 
PITCH(N°K)=Y C 1 ) 

IF(NACT.EO.O) GOTO 16 
on io tac=i#nact 
FCIAC,NPW)sO.O 
DO 15 J*=1,ISS 

F(IAC,NPK)=F(I ACfNPK)+AK(lAC, JK)*Y{ JK1 

CONTINUE 
00 19 JKNsWISS 

F(IAC#NPK)sF(IAC f NPK)+AKf IAC/ JKMf ISS) *Y(JKM+NMODES) 

) CONTINUE 

4 CONTINUE • ... 

> CONTINUE 

00 0 J=1,NP0NTS 
Of JlsO.O • 

00 3 I =2 » NMODES 

COSH.^s (FXP(FREO( I)) +EXP(-FREQ(I)n/2. 

SINH,v = (cxP(FREOtin-EXP(-FREOCI)n/2. * * 

C0SH'NL=(EXP(FREG(n*AL(Jn*EXPC-FPEQ(I)*AlCJn)/2. 
SINHAL=CEXPCFREQ(n*AL(Jn«*EXP(-FREQ(n*ALf J) ))/2. 

OUPl = (COSCFREO(in-COSHW)/(SINHW-SIN(FREam n 

DUP2=SINCFRE0(n*ALCJ)n-SINHv«L 

DUP3=C0S(FREGfI ) *AL ( Jl) ♦COSHWU 

0UP=0UP1*0UP2‘DUP3 

0(J)=0( J)fYCI) *DUP 

CONTINUE 

continue 

0MAx=gn ) 

OMiMsQU ) 

00 5 Msl,MPOMTS 
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IF(Qmax.LT.QCM+1) ) goto 6 • * 

GOTO 7 ** 

> QWAXsOfWl ) 

r IF(QMIN.LT.Q( w +-t ) ) GOTO 5 

QMINsQ f M+ 1 ) 

» CONTINUE 

WRITEC8, 18) Qmax.QMIM 

8 FORMA TCPOX, '0 M AX=',E13.6, 'QMIN=',E13.6, '**iw THIS CASE') 

IF(ABSCQmax) .GT.AOSCQMIN)) GOTO 8 
QMAXsAHSCOVIN) 

GOTO 9 

1 QMIN=-ABS(Qmax) 

' CONTINUE 

IFCABSfOMAX) .LT.OSMALL) GOTO 12 
IF(QMAX.Gr.O M AXl) GOTO 1910 

QMAXzQMAXl 
QMINs-qmax 1 
GOTO 1911 

910 QMAX1=QMAX 

911 CONTINUE 
WQITE(8,I0) QMAX,0HIN 

0 FORMAT Cl OX, '**G)MAX=',E13.&,5X, '**0MIN=» ,E13.6) 
IF(NSELCT.EO.O) GOTO 17 
CALL PLOGOCO.O, 1 .25) 

CALL PSIZEC8. 0,1.0) 

CALL PG 5 1 0(0,2) 

7 CONTINUE 

CALL PLOX (QMIN,Q,QVAX,NPnNTS) 

2 CONTINUE 
CONTINUE 
RETURN 
END 
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